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■ We theoretically study the effects of elastic collisions on the determination of frequency standards 

<^^> via Ramsey fringe spectroscopy in optical-lattice atom clocks. Interparticle interactions of bosonic 

atoms in multiply-occupied lattice sites can cause a linear frequency shift, as well as generate 
asymmetric Ramsey fringe patterns and reduce fringe visibility due to interparticle entanglement. 
' We propose a method of reducing these collisional effects in an optical lattice by introducing a phase 

difference of n between the Ramsey driving fields in adjacent sites. This configuration suppresses site 
to site hopping due to interference of two tunneling pathways, without degrading fringe visibility. 
qq ' Consequently, the probability of double occupancy is reduced, leading to cancellation of collisional 

i ' shifts. 
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PACS numbers: 42.50.Gy, 39.30.+W, 42.62.Eh, 42.50.Fx 



I. INTRODUCTION 



Current state-of-the-art atomic clock technology is based mainly on trapped single ions or on clouds of free falling 
cold atoms [j. Recently, a new type of atomic clock based on neutral atoms trapped in a deep "magic- wavelength" 
c/2 ■ optical lattice, wherein the ground and excited clock states have equal light shifts, has been suggested 0,0- For 
interfering, far-detuned light fields, giving rise to a "washboard" intensity pattern, the optical potential experienced 
by an atom in a given internal state is Vo(r) = H\Qo(r)\ 2 / (4A ) where £lo( r ) = 2^i? (r)/ft is the dipole coupling 
frequency, Eo(r) is the radiation electric field strength at position r, and Ao is the detuning of the oscillating optical 
field from resonance. In a magic-wavelength optical lattice, the same exact potential is experienced by the ground 
and the excited state of the clock transition. 

Optical-lattice atomic clocks are advantageous due to the suppression of Dopplcr shifts by freezing of translational 
motion (the clock operates in the Lamb-Dickc regime, where atoms are restricted to a length-scale smaller than 
the transition wavelength Q), the narrow linewidth due to the long lifetime of the states involved in the clock 
transition, and the large number M of occupied sites, minimizing the Allan standard deviation Q]. As already 
mentioned, the optical-lattice potential light shift is overcome by using trapping lasers at the magic wavelength 0, U . 
, Operating optical-lattice clocks in the optical frequency regime, rather than in the microwave, has the added benefit 
of increasing the clock frequency v and thereby reducing 5u/u by five orders of magnitude. A coherent link between 
optical frequencies and the microwave cesium standard is provided by the frequency comb method Q . Thanks to these 
■ characteristics, atomic optical-lattice clocks have promise of being the frequency standards of the future. Attempts 
C/3 | to further improve the accuracy of such clocks using electromagnetically induced transparency to obtain accuracies 
on the order of 10 17 or better have also been suggested 0. 

A 3-dimensional (3D) optical lattice configuration would allow the largest number of atoms that can participate in 
an optical-lattice atom clock. However, when sites begin to be multiply-occupied, atom-atom interactions can shift 
the clock transition frequency This is particularly problematic in a very deep optical lattice since the effective 
density in sites with more than one atom will then be very high due to the highly restricted volume. Hence, the 
collisional shift, proportional to the particle density, can be very large. It is therefore important to ensure that atoms 
individually occupy lattice sites (preferably in the ground motional state of the optical lattice). One way of achieving 
this goal is by low filling of the optical-lattice, so that the probability of having more than one atom per site is small. 
If hopping of an atom into a filled lattice site during the operation of the clock is small, collisional effects will be 
minimal. 

One kind of optical lattice clock configuration that can have minimum collisional effects uses ultracold spin-polarized 
fermions in a deep optical lattice. For example, a system of ultracold optically pumped Sr atoms in the 5s 2 1 So 
\F = 9/2, Mp = 9/2) internal state, filled to unit filling. The nuclear spin of the atoms is protected and one expects 
that the gas will remain spin polarized for a very long time. The probability of finding two fermionic atoms in a deep 
lattice site is extremely low in such a system as long as the gas is sufficiently cold, since higher occupied bands will 
not be populated, and only one fermionic atom will be present in a given site due to the Pauli exclusion principle. 
Moreover, fermionic atoms cannot interact via s-wave collisions and p-wave and higher partial waves are frozen out 
at low collision energies. Provided the ground motional state of the spin polarized fermion system can be attained by 
adiabatically turning-on the optical lattice |lfj | , this system seems to be extremely well-suited for making an accurate 
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clock. 



Another potentially interesting configuration involves ultracold bosonic atoms such as 88 Sr atoms in the 5s 2 1 So 
\F = 0,Mp = 0) internal state in an optical lattice with very low filling. Since the atoms are bosons, there is no 
single- atom-occupancy constraint, allowing for collisional shifts. These shifts can be minimized if the filling factor p is 
small. Moreover, multiple occupancy caused by tunneling between adjacent populated sites is reduced if the optical 
lattice is deep and the probability of hopping of an atom into an adjacent filled site, Jhop, is small. Therefore, we 
expect that the collisional shift should be very low. Here, we investigate the Ramsey fringe clock dynamics for such a 
system. The transition between the 1 So \F = 0, Mp = 0) and the 3 Po \F = 0, Mf = 0) state of 88 Sr can be a Raman 
transition but we can think of the transition between 1 So and 3 Po as a Rabi transition as long as the detuning 
from the intermediate state coupling these two levels is large. Yet another possibility is to use a weak static magnetic 
field to enable direct optical excitation of forbidden electric-dipole transitions that are otherwise prohibitively weak 
by mixing the 3 Pi with the 3 Po state @. In contrast to multiphoton methods proposed for the even isotopes 
this method of direct excitation requires only a single probe laser and no multiphoton frequency-mixing schemes, so 
it can be readily implemented in existing lattice clock experiments Q. I.e., this method for using the "metrologically 
superior" even isotope can be easily implemented. However, one of the problems that can arise is that more than one 
bosonic atom can fill a lattice site, and these atoms can interact strongly. We shall assume that sites are initially 
populated with at most one atom per site, but that during the operation of the Ramsey separated field cycle, i.e., the 
delay time between the two 7r/2 Ramsey pulses, T, atoms from adjacent sites can hop; if a site ends up with more 
than one atom during the cycle, the clock frequency will be adversely affected by the collisional shift. 

It is easy to obtain an order of magnitude estimate for the collisional shift 8v in this kind of clock. The product of 
the hopping rate Jh op and the Ramsey delay time T gives the probability for hopping between sites. As will be shown 
below, the shift obtained is therefore given by 5v = J\ W pTQ, where Q is the nonlinear interaction strength parameter 
defined in Eq. I|18|) . If we consider only a single site, then to attain an accuracy of one part in 10 17 , one must have 
Jho V TG < 10~ 17 v, whereas for N s occupied sites, AopTQ / y/N~ s < KT 17 ia The precise values of these parameters can 
vary greatly for different experimental realizations. In particular, the hopping rate Jhop is exponentially dependent 
on the optical lattice barrier height and on particle mass. 

Three distinct collisional effects are found for a single, multiply-populated lattice site. These include a simple linear 
frequency shift, a nonlinear shift resulting in asymmetric fringe patterns, and an entanglement-induced reduction in 
fringe contrast. In order to show how these effects can take place dynamically during the application of the Ramsey 
pulse sequence, we consider a ID optical lattice and focus on two adjacent sites filled with one atom in each site, 
and calculate the probability of double occupancy due to hopping (tunneling) of an atom to its adjacent site and 
the resulting collisional shifts. We show that by varying the direction of the Rabi drive laser with respect to the 
principal optical lattice axis, one can induce a phase-difference between the Rabi drive fields in adjacent lattice site. 
We find that due to interference of tunneling pathways, hopping is suppressed when this phase difference is set to n, as 
compared to the case where all sites are driven in phase. Consequently, detrimental collisional effects arc diminished 
for inverted-phase driving of adjacent sites. 

The outline of the paper is a follows. In Sec. |H] we review the Ramsey separated oscillatory fields method for 
noninteracting particles and set up the notation we use in what follows. In Sec. IIIII we construct a second-quantized 
model, describing the single-site Ramsey scheme for interacting bosons. This model is numerically solved in Sec. II VI 
demonstrating various collisional effects on Ramsey fringe patterns. Hopping between sites is introduced in Sec. 
where we present two-site results and propose an inverted Rabi phase scheme to cancel collisional shifts. Conclusions 
are presented in Sec. IVII 



Norman Ramsey introduced the method of separated oscillatory fields in 1950 |ll|. A long time-period between the 
application of two nearly resonant coherent fields makes the Ramsey resonance very narrow, and thus suitable for high- 
performance clocks and precision measurements |l2l 1 1 3| . The method has since become a widely used technique for 
determi ning resonance frequencies to high precision. For example, in the Cs fountain clock experiments summarized 
in Refs. \l3L Il4j . the observed linewidth of the Ramsey resonance was 1 Hz, two orders of magnitude below that of 
thermal Cs clocks H3, EE 

For a two-level atom in an intense short near- resonant pulse with central frequency to, the Hamiltonian in the inter- 
action representation with \ip) = a g (t) exp(— i(e g / 'h+ uj)t)\g) +a e (t) exp(— ie e t/h)\e) , and rotating wave approximation, 
takes the form 



II. RAMSEY SEPARATED OSCILLATORY FIELDS 




(1) 
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FIG. 1: (color online) Collisional shift and loss of single-particle coherence in a Ramsey scheme. The solid red line traces the 
trajectory of the Bloch vector in the absence of interactions. When interactions are present (solid blue line), the frequency of 
phase oscillations is modified if G gg 7^ G ee and the length of the Bloch vector is not conserved when G gg + G ee 7^ 2G ge (note 
that the solid blue curve deviates from the dashed great circle during the free induction decay, whereas the solid red curve does 
not). 



where ft = 2fj,A/h is the Rabi frequency, A is the slowly varying envelope of the electric field strength, [i is the 
transition dipole moment, and A = (e e — e g — huj)/H is the detuning from resonance of the laser frequency lo. The 
solution of the optical Bloch equations for the two-level atom is given in terms of the unitary evolution operator for 
the two-level system for a real slowly varying envelope: 

/ ™s(%V f^in(^) gsin(%^) \ 

U = . (2) 

V gsm(¥) cos(^) + fA si n(%^ 



4 




Here fi s = -\/|f2| 2 + A 2 is the generalized Rabi frequency. In the Ramsey method, the system, assumed to be initially 
in the ground state \g), is subjected to two pulses separated by a delay time T, 

( n if o < t < t p , 

fi(£) = \ if t p < t < T + T p , (3) 
{ n if T + Tp < t < T + 2t p , 

with fir p = it/2 and T 3> t p . From transformation J5J), it is clear that the effect of the first pulse is to evolve the 
initial ground state \g) into the state (\g) + i\e))/y/2. In a Bloch-sphere picture with u = di(a*a e ), v = 3(a*a e ), and 
w = (\a e \ 2 — \a g \ 2 )/2, the Bloch vector (u, v, w) is projected by the first pulse into the uv plane, as depicted by the 
red line in Fig. [J] During the delay time, the system carries out phase oscillations, corresponding to rotation of the 
Bloch vector in the uv plane with frequency A. Finally, the second pulse rotates the vector again by an angle of Q g T p 
about the u axis, as shown in Fig.^ Fixing A and measuring the final projection of the Bloch vector on the w axis as 
a function of T, one obtains fringes of fixed amplitude Q/Q g and frequency A. Alternatively, fixing T and measuring 
w(t > T + 2t p ) as a function of the detuning A, results in a power-broadened fringe pattern of amplitude f2/f2 g and 
frequency 2-k/T. The resulting probability to be in the excited state is given by 

P - = 7ST sin2 (^9 T p) (~ C0 <ln g T p ) cos(iAT) + A sm(ln g T p ) sin^AT)) . (4) 



n 2 v 2 flP'V 2 2 "a 2 2 / 

Figure |5] shows the population inversion w versus time and detuning A using a Ramsey separated field method 
for one atom in an optical lattice site. The final time corresponds to the time at which the Ramsey clock signal 
is measured as a function of detuning A, i.e., either the population of the ground or the excited state is measured 
as a function of A. Note that the excited state population at the final time is unity at zero detuning and that the 
population inversion oscillates as a function of detuning. 

It is easy to generalize this treatment to a time-dependent Rabi frequency U(t) = 2fiA(t)/H due to a pulse of light 
which turns on and off with a finite rate in terms of the pulse area, JJ" P Q(t')dt' . For a tt/2 pulse, JJ" P Q,(t')dt' = tt/2, 

where r p is the pulse duration. For the Ramsey pulses, J Q P tt(t')dt' = tt/2, J t v Q(t')dt' = 0, and J T+T p Q,(t')dt' = 
tt/2. 



III. SECOND-QUANTIZED RAMSEY MODEL 

In order to study the effects of collisions on Ramsey fringes obtained in a separated-fields scheme, we use a second- 
quantized formalism to treat a multiply-populated single site of the optical potential and calculate the Ramsey fringe 
dynamics. The many- body Hamiltonian for a system of two-level atoms, all in the same external motional state of a 
trap, can be written as [2(J, l21j 



m(t) 

2 { 

t=9,e i,3=9,e 



H=^2,E i a\a i — {a\a e + a\a g ) + ^ G tj d\d]a :j d l . (5) 
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Here the subscripts g and e indicate the ground and excited states of the two level atom, and hi is the annihilation 
operator for an atom in internal state i, where i = g,e. The self- and cross- atom-atom interaction strengths are 
denoted as Gu and Gij respectively, the internal energy of state i is denoted as Ei, where E g — e g + hw and E e = e e . 
For bosonic atoms, the creation and annihilation operators satisfy the commutation relations 

[a u at] = % , [al,a]} = , [a,, aj] = . (6) 

Defining the operators 



~27 



ni = aja,-, iV = ?i g + n e , L 2 = ^— , = , L y = — — , (7) 



the Hamiltonian can be written as 

2 

H = EN + HAL Z - m(t)L x + G'ti {n\ - n q ) + 2G ge n g n e , (8) 

t=i 

where 

E = (fi^ + e 3 + e e )/2, (9) 
HA = {hw-€ e + e g ). (10) 

Since N commutes with H, the Hamiltonian conserves the total number of particles. Typically, in a Ramsey- fringe 
experiment, the initial state is assumed to have all N atoms in the ground state \g). Fixing N = N in the single 
site that we are considering here (i.e., so that there are no fluctuations), and using the identities n g = N /2 + L z and 
h e = N/2 — L z , we finally obtain 

C 4- C M 2 
H = {E _^ l +^ )N + {Ggg + Gee + 2Gge) ^_ 

-m(t)L x + [HA + (G gg - G ee )(N - 1)}L Z 

+ {G gg + G ee - 2G ge )L 2 z . (11) 

Since Ramsey spectroscopy measures essentially single-particle coherence, we will be interested in the dynamics of 
the reduced single-particle density matrix (SPDM) p{t): 



N 

p(t) = (p{t)) = —2 + u(t)a x + v(t)a y + w{t)a z , (12) 



where p(t) = a](t)dj(t) is the reduced single-particle density operator, 2 is a two-by-two unity matrix, crj, i = x, y, z are 
Pauli matrices, and u = (L x ), v = (L y ), w = (L z ) are the components of the single-particle Bloch vector, corresponding 
to the real- and imaginary parts of the single-particle coherence, and to the population imbalance between the two 
modes, respectively The expectation values, (•) are over the iV-particle states. The Liouville-von-Neuman equation 
for the SPDM, 

±p=l([p,H}) = ^ Yl "J" . ( 13 ) 

j=x,y,z 

is thus equivalent to the expectation values of the three coupled Heiscnberg equations of motion for the SU(2) 
generators, 

jl x = — [A + U(N - l)]L y - Q(L y L z + L z L y ), (14) 
jl y = [A + U{N -l)]L x + Q{L x L z + L z L x ) 

+n(t)L z , (15) 

^L z = -n(t)L y , (16) 
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FIG. 3: (color online) Linear collisional shift. Ramsey fringes for N particles in a lattice site, for (a) U = Q = and (b) 
U(N - 1) = 3000 Hz, g = 0. 



where we denote the linear and nonlinear interaction strength parameters respectively by 

U = (G gg - G ee )/h , 



(17) 



Q = (G 



99 



Gp. 



2G ge )/h . 



(18) 



In order to study effects resulting in from the coupling to an environment consisting of a bath of external degrees 
of freedom, we use the master equation [lgllisj 



dt' 



{o{6 k ,p}), 



(19) 



where the second term on the right hand side of l|19|) has the Markovian Lindblad form ^7J an d gives rise to dissipation 
effects of the bath on the density operator. Such terms may be used to depict decay due to spontaneous emission, atom- 
surface interactions, motional effects, black-body radiation and other environmental effects. The Lindblad operators 
Ok are determined from the nature of the system-bath coupling and the coefficients T k are the corresponding coupling 
parameters. In what follows we shall assume that the most significant dissipation process is dephasing, and we neglect 
all other dissipation effects (e.g., we assume spontaneous emission is negligible because the lifetime of the excited 
state is much longer than the Ramsey process run-time, etc.). Thus, the Lindblad operator is taken to be L z , yielding 
dephasing of the transition dipole moment without affecting the population of the ground or excited states. We shall 
only study this type of ti dephasing, without fully exploring the effects of other types of Lindblad operators. Clearly, 
it is easy to generalize this and consider the effects of additional Lindblad operators, but we do not do so here. 
Moreover, non-Markov treatments can be used to generalize the Markovian approximation made in deriving Eq. 11911 
111 111. 
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IV. SINGLE-SITE MANY-BODY DYNAMICS 

We first study how interparticle interactions can modify the Bloch-vector dynamics in a single-site Ramsey scheme. 
From Eq. and Eqs. JTU'dU it is evident that there are three possible collisional effects. First, differences between 
the interaction strengths of ground and excited atoms, will induce a linear frequency shift of the center of the Ramsey 
fringe pattern by an amount U(N — 1). This effect is illustrated in Fig.^ where the trajectory of the Bloch vector 
during a Ramsey sequence is traced in the absence of interactions (red) and in their presence (blue). Clearly, the 
frequency of oscillation in the uv plane during the delay time T is modified by the interaction. In Fig. |21we plot the 
resulting fringe pattern when Q is set equal to zero. Fig. Ofa) depicts the Ramsey fringes without any interaction, 
hi = Q = 0, whereas in Fig. E^b) we plot the fringe pattern with hi = 3000 Hz. Since for Q = Eqs. (|14fl - l|16[l remain 
linear, and the only effect is a linear frequency shift. The whole curve w(A) versus A is simply shifted in A by an 
amount U{N - 1). 

Also shown in Fig. Q] is the loss of single-particle coherence due to interparticle entanglement [2(j ■ This dephasing 
process is only possible for the nonlinear case when Q ^ 0. It is manifested in the reduction of the single-particle 
purity p = Tr(p 2 ) = \u\ 2 + \v\ 2 + \w\ 2 , during the evolution. In comparison, for the single-particle case depicted in red, 
single-particle purity is trivially conserved and the length of the Bloch vector is unity throughout the propagation 
time; the Bloch vector at the end of the process sits well within the Bloch sphere. 

The decay of single particle coherence due to entanglement is illustrated in Fig. where the Bloch vector (u, v, 
w,), the single particle purity p, the single-particle entropy s = Tr(p\np) = — ln[(l + v / p) 1+ v / P(l — v ^p) 1_ v^/4]/2, 
and the variance of w, Aw, are plotted versus time in a Ramsey separated field method for A = 0. The total time 
for the Ramsey process is taken to be t tot = T + 2t p = 1.0 x 10 -3 s, with the first and second Ramsey pulses each 
of duration t t ot/12, and the Rabi frequency of the pulses are 3.0 x 10 3 Hz, so the two pulses each have pulse area 
7r/2. The interaction strengths were arbitrarily set to G ee = G ge = 0, so that hi = Q = G gg /h = n. In frame (a) k is 
set to zero (corresponding to the case of one atom per site in the optical lattice) . The effect of decoherence due to 
the Yl z term in Eq. (|13|l is shown in frame (b), while collisional dephasing is depicted in frames (c)-(e) were we set 
k = 100 Hz. In frames (c) and (e) we have taken two atoms per site, whereas in frames (d) and (f) we have taken 10 
atoms per site. In frames (e) and (f), dephasing with strength 7 Z = 1.0 x 10 2 Hz is included. From comparison of 
Fig- Elk) and Fig.^c-d), it is clear that the loss of single-particle coherence due to entanglement is similar to the effect 
of dephasing due to coupling to an external bath. In particular, the final population inversion is strongly affected. 
Comparing Figs.^Jb) and Fig.^Jc), it is clear that collisional dephasing is stronger as the number of particles grows. 
The combined effect of decoherence and interparticle entanglement is shown in frames (e) and (f). 

Yet another collisional effect is caused by the QL 2 term in the Hamiltonian l|ll|) . In the mean-field approximation, 
replacing Li by their expectation values, this term leads to a nonlinear frequency shift of 2Qw{t p ), where w(t p ) is 
the value of the population imbalance after the first pulse. As illustrated in Figs. [3] and the nonlinear shift results 
in an asymmetric fringe pattern, because for G ^ 0, w(t p ) depends on the sign of A during the first pulse. Without 
interactions, when li = Q = 0, the dynamical equations (|14|> - i|16fl are symmetric under A — > —A , L x — > — L x . 
This symmetry shows up in Fig. |3{a) where the evolution with two opposite values of A is traced over the Bloch 
sphere. As shown in the various projections onto the uv, vw, and uw planes, ma(^) = — u -a(0> ^aW = ^-aW, an d 
WA(t) = W-a(^), where (u^(t\v^(t\w^(€)) denote the Bloch vector (u(t),v(t),w(t)), evolved in a Ramsey sequence 
with detuning A. Similarly, when Q = and hi is nonvanishing, the resulting shifted pattern is still symmetric about 
A = —(A — 1)U since the dynamical equations are invariant under A + (N — 1)U — > —(A + (N — 1)U) , L x — > — L x . 
However, for finite values of Q the symmetry breaks down because the last terms on the r.h.s. of Eqs. Ijl4|l an d 
(II 51) change sign, as demonstrated in Fig.|S{b). Consequently, when the final projection on the w axis is plotted as a 
function of the detuning A, we obtain asymmetric fringe patterns. In Fig. [HJ we compare the symmetric single-particle 
pattern (Fig. |SJi) with asymmetric two-particle Ramsey fringes obtained with nonvanishing values of Q (Figs. 03 and 
El;). The linear shift U is set equal to zero for all plots. Both the asymmetry of the emerging lineshape and the 
reduction in fringe contrast due to collisional dephasing are evident. 

The combined effect of the linear shift and the nonlinear term in Eqs. H14f> - I|16|) i s shown in Fig. |7| where the final 
population inversion w(t = r + 2r p ) is plotted as a function of the detuning A and the interaction strength, set 
arbitrarily to hi = Q = k. Fringe maxima arc located about (A + k(N — 1) + 2kw(t p )) = 2jn where j is an integer and 
w(t p ) depends on both A and k. As k increases, collisions shift the peak intensity and the fringe contrast is reduced 
due to collisional dephasing. 
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FIG. 4: (color online) Bloch vector, u, v, w, purity p, entropy s and variance of w versus time using a Ramsey separated 
field method for zero detuning A. (a) Two atoms without dephasing and without interaction, (b) Two atoms with dephasing 
(r_L. = 100 Hz) and without interaction, (c) Two atoms without dephasing but with interaction (k = 100 Hz), (d) Ten atoms 
without dephasing but with interaction (k — 100 Hz), (e) Two atoms with dephasing (Tl z = 100 Hz) and interaction (k = 100 
Hz), (e) Ten atoms with dephasing and interaction. 



V. TWO-SITE MANY-BODY DYNAMICS 



Having established how intcrparticlc interactions affect the Ramsey lineshapes in a multiply-occupied single lattice 
site, we proceed to consider an optical lattice with multiple sites labeled by the site index a. Each site is populated 
with atoms that can be in any one of two levels, the ground state level \g) and the excited state level |e), coupled 
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FIG. 5: (color online) Collisional asymmetry of Ramsey fringes. Red and blue curves correspond to Bloch-vector trajectories 
with positive and negative detuning, respectively. The final position of the Bloch vectors is indicated by a dot at the end of the 
trajectories. In the absence of interactions (a) the final projection onto the w axis is independent of the sign of A. However, 
when interactions are present (b) nonlinear collisional shifts break the symmetry. Dash-dotted lines depict projections of the 
trajectory onto the uv, uw, and vw planes. 

through a Rabi flopping term in the Hamiltonian. The full Hamiltonian is given by 

H = E x&f ( 20a ) 



rr(a) , fr( a ) , fr(a) , fr( a ) 
-"int "T" -"diff -"sum T "R 



^hop^ = Who P [(4«/3 + h.c.) + {btfip + h.c.)} , (20b) 
H$ = ^U aa al l al 1 a a a a + ^U bb P a P a b a b a 

+Uabal c bl l b a a a , (20c) 



HA 



m = -Y(ala a -bib a ) , (20d) 



- (a + hw + 6b fat a +tfb ) 

-"sum r, V"q u Q i u a u a) j 



(20e) 



2 

4 a) = ^(aLSae^-W+Staae-^-W). (20f) 

Here, the operators a and b are bosonic destruction operators for atoms in the two states \g) and |e), and the indices 
a and (3 run over the lattice sites, with (a, /?) denoting a pair of adjacent sites. The constants Jhop and E/y , i, j = g,e, 
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FIG. 6: (color online) Asymmetric Ramsey fringes. Numerically calculated Ramsey fringe patterns, at three different values of 
the nonlinear interaction: (a) Q = 0, U = 0, (b) Q = 1500Hz, U = 0, and (c) Q = 3000 Hz, U = 0. 
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FIG. 7: (color online) Population inversion w, versus detuning A and interaction strength U — Q — « in a Ramsey separated 
field configuration at the final time after the second pulse. 
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FIG. 8: (color online) Probability for double occupancy P2, during a Ramsey sequence with A = 0, versus relative Rabi-drive 
phase angle 6 for Jh op = 100 Hz. Tunneling is suppressed for 9 = n due to destructive interference and enhanced for 6 = due 
to constructive interference. Interparticle interactions are set equal to zero. 



are the strength of the hopping to neighboring sites in -ffh op and the effective on-site interaction in H lnt , respectively, 
-ffdiff and -ffsum are the energy difference and average of the dressed states of the atoms, and Hr induces Rabi 
transitions between the two atomic internal states with Rabi frequency Q.{t) and phase <fr a (t), which are related to 
the intensity and phase of the laser which induces these transitions. The Hamiltonian in i|20|) is identical to what we 
used in previous sections, except for the addition of the hopping term -ffhop that can result in hopping of particles to 
adjacent sites with rate Jhop- 

Note that the phase 4> a (t) can depend upon the site index a. Consider, for simplicity, a ID optical lattice along 
the x axis, and a plane wave field with detuning A from resonance with wave vector k. The Rabi frequency at site 
a is given by fJ(r Q ,£) = f2(t)exp[ik • r a ] = Q,(t) exp(ik x x a ). Thus, there is a phase difference 9 between the Rabi 
frequency at neighboring sites, 9 = k x Sx, where Sx is the lattice spacing. The angle between the wave vector k and 
the x axis can be adjusted to control the phase difference 9 — (f> a — 4> a -i- In what follows, we show that a proper 
choice of the relative phase angle 9 between Rabi drive fields in adjacent sites may be used to suppress the tunneling 
between them and thus reduce collisional shifts. When Jh op = 0, there is no relevance to the phases (f> a {t) (they do 
not affect the dyanmics), but as soon as hopping from one site to an adjacent site is allowed, these phases can affect 
the dynamics. 

We expect that for sufficiently low densities, the probability of finding more than two adjacent populated sites is 
very small. It is therefore possible to capture much of the physics of the hopping process, using a two-site model. 
At higher densities, more elaborate models should be used. In Fig. ISlwe plot numerical results for two particles in 
two sites (N = 2, a = 1,2) for Ji lop = 100 Hz, showing the probability of double occupancy during a Ramsey fringe 
sequence in the presence of tunneling, as a function of the relative angle 9. The parameters used are the same as 
those used previously, except that now we allow for hopping, and therefore Jh op and the angle 9 affects the dynamics. 
The initial conditions correspond to a single atom in its ground state in each lattice site before the Ramsey process 
begins. When both atoms are driven in-phase (9 = 0), tunneling takes place, leading to multiple occupancy. The 
tunneling is significantly suppressed for a 9 = 7r phase difference between the Rabi drives in adjacent sites. As shown in 
Fig- El this suppression originates from destructive interference. In Fig. [3]} we plot the level scheme for a two-particle, 
two-site system. The levels |nf , n\ ; nf, n|) denote Fock states with n 9 a and n e a particles in the ground- and excited 
states respectively, of site a. Depending on the relative phase between the optical drive fields, tunneling from states 
with a ground state atom in one site and an excited atom in another, can interfere constructively or destructively. 
Therefore, as seen in Fig. Eb and Fig. Et, tunneling can only take place between even-parity states. For in-phase 
Rabi drive 9 = (Fig. Eb)i the initial even-parity state |1,0; 1,0) is coupled by a single optical excitation to the 
even-parity state (l/v2) (|0, 1; 1, 0) + |1, 0; 0, 1)), for which the two tunneling pathways interfere constructively. In 
contrast, for out-of-phase drive 9 = tt (Fig.Et), the first tt/2 pulse drives the system partially into the odd-parity state 
(l/\/2) (|0, 1; 1, 0) — |1, 0; 0, 1)) which is 'dark' to tunneling, leading to a reduced probability of multiple occupancy 
during the time evolution. 
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I2,0;0,0> ll,0;l,0> I0,0;2,0> 



I0,2;0,0> - I0,0;0,2> I0,1;0,1> I0,2;0,0> + I0,0;0,2> 




I0,0;2,0> - I2,0;0,0> ll,0;l,0> I0,0;2,0> + I2,0;0,0> 



( c ) I0,2;0,0> - 10,0;0,2> I0,1;0,1> I0,2;0,0> + 10,0;0,2> 




I0,0;2,0> - I2,0;0,0> ll,0;l,0> I0,0;2,0> + I2,0;0,0> 



FIG. 9: (color online) Two-particle level schemes for Ramsey spectroscopy in the presence of tunneling: (a) general scheme (b) 
coupling between parity eigenstates for 9 = (c) coupling between parity eigenstates for 8 = n. Rabi coupling between ground- 
and excited states, is denoted by red and green solid arrows, corresponding to different phases of the driving fields in adjacent 
sites. Blue dashed arrows denote hopping between sites. Fock states are denoted as nf , nf; nf , nl), wh ere n 9 a = (aj,a a ) and 
< = {bt,b a ). 
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FIG. 10: (color online) Probability of double-occupation at the end of a Ramsey sequence versus detuning, for different values 
of the nonlinearity U = Q = k: (a) k — 0, (b) n = 5.0 x 10 2 , and (c) k = 1.0 x 10 3 . Red curves correspond to 9 = whereas 
blue curves show the lineshape for 6 = n. 



The lower probability of finding multiply occupied sites in an inverted-phase Rabi-drive configuration, is manifested 
in the emerging Ramsey fringe patterns. In Fig. llOl we plot the probability of double-occupancy at the end of a Ramsey 
sequence, as a function of detuning, for different degrees of nonlinearity which is arbitrarily set to U = Q = n. While 
increasing nonlinearity generally tends to localize population and reduce tunneling, the probability of finding both 
particles in the same lattice site is always lower for out-of-phase driving. The resulting Ramsey lineshapes are plotted 
in Fig. 1111 clearly demonstrating reduced collisional effects for 8 = it (blue curves), compared to the in-phase driving 
fields scheme (red curves). The center of the distributions for 9 — are clearly strongly shifted from A = 0, and 
much less shifted for 9 = tt. This reduced shift in the center of the distribution is, of course, of central importance 
for a clock. 



VI. CONCLUSION 



Optical lattice atom clocks hold great promise for setting frequency standards. In order to achieve high accuracy 
with boson atoms, collisional frequency shifts have to be taken in account. We have shown that collisions can degrade 
Ramsey lineshapes by shifting their centers, rendering them asymmetric, and by reducing fringe-visibility due to the 
loss of single-particle coherence. Considering two adjacent populated sites in a ID optical lattice configuration, we 
propose a method of reducing dynamical multiple population of lattice sites, based on driving different sites with 
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FIG. 11: (color online) Ramsey fringe patterns for different values of the nonlinearity U — Q — k: (a)« = 0, (b) k = 5.0 x 10 2 , 
and (c) k = 1.0 x 10 3 . Red curves correspond to 8 = whereas blue curves show the lineshape for 9 — ty. The arrows serve to 
roughly indicate the centers of the distributions. 



different phases of the Rabi drive. Due to destructive interference between tunneling pathways leading to states with 
one ground and one excited atom in the same site, hopping is reduced and collisional effects arc canceled out. It is 
difficult to see how to implement this kind of interference in a 2D or 3D configuration. 

While collisional effects are considered here as an unwanted degrading factor, our work also suggests the potential 
application of the Ramsey separated field method for studying entanglement. 
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